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This review gives an overview of most recent measurements of top quark production cross sections 
including differential cross sections and searches for new physics in the top quark sector. Datasets 
corresponding to an integrated luminosity of up to 5.7 fb~^ are presented which were taken at the 
Tevatron proton-antiproton collider at a center-of-mass energy of ^/s = f .96 TeV at Fermilab. 

O ■ I. TOP QUARK PRODUCTION AT THE TEVATRON 

' Top quarks were first observed via pair production at the Fermilab Tevatron Collider in 1995 The Tevatron 
is still the only place where top quarks can be studied with a very high precision. At the Tevatron, top quarks arc 
either produced in pairs via the strong interaction or singly via the electroweak interaction i2|]. In the framework 
of the standard model (SM), each top quark is expected to decay nearly 100% of the times into a W boson and a 
^ ' ' h quark Q . W bosons can decay hadronically into qq' pairs or leptonically into eve , nv^j. and tvt with the r in turn 
^ ' decaying into an electron, a muon, or hadrons, and associated neutrinos. 

■ In ti production, if both W bosons decay hadronically the final state is called all-hadronic (or all-jets) channel. 
O If one of the W bosons decays hadronically while the other one produces a direct electron or muon or a secondary 



> 



electron or muon from t decay, the final state is referred to as the .^-|-jets channel. If both W bosons decay leptonically, 
this leads to a dilepton [U) final state containing a pair of electrons, a pair of muons, or an electron and a muon, or 
a hadronically decaying tau accompanied either by an electron or a muon (the r£ channel) . 



■ II. TOP QUARK PAIR PRODUCTION 

(N 



Exploring the top cross section in different decay channels and using different assumptions is important because 



■ signs of new physics might appear differently in the various channels. 

• i— ( , A. Top Quark Pair Production Cross Section 

>< 

' The top quark pair production cross section cr^j is known to high accuracy in the SM The measurement 

of the top quark pair production cross section therefore provides an important test of calculations in higher order 
Quantumchromodynamics (QCD) including soft gluon resummations. Any deviation from the SM prediction of the 
measured ti cross section could either be a hint for new physics in top quark pair production or in top quark decays. 
For example, an exotic decay of the top quark, such as the decay into a charged Higgs boson and a h quark [t — > H^h) 
would lead to deviations of the measured in individual final states compared to the SM prediction. 

Usually there are two different techniques to enhance the top pair signal over the background. The "topological" 
method explores tt event kinematics to distinguish tt signal from background. Because of the large mass of the top 
quark, ti events are more energetic, central, and isotropic compared with the dominant backgrounds such as W^+jets 
and QCD multijet events, whose kinematics are more influenced by the boost from the momentum distribution of 
the colliding partons. The "counting" method utilizes the identification of jets originated from 6 quarks (6-tagging). 
The first method does not rely on the assumption that a top quark decays into a 6-quark as opposed to the second 
one. Thus they are sensitive to different systematic uncertainties. 
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One example for a "topological" analysis from the DO Collaboration [9[ uses an integrated luminosity of 4.3 fb~ 
exploring final states with 2, 3 or > 4 jets, thereby defining twelve disjoint data sets. To distinguish tt signal from 
background, a discriminant is constructed that exploits differences between kinematic properties of the tt f+jets 
signal and the dominant M^+jets background. A multivariant discriminant function is calculated by a random 
forest (RF) of boosted decision trees (BDT). The output of the BDT discriminant for > 4 jets is presented in 
Fig. [T] (left). To extract the tt cross section one performs a binned maximum likelihood fit of the discriminant 




Discriminant Output 



FIG. 1: Left: output of the BDT discriminant for £+> 4jets, showing contributions from tt signal using a cross section of 
7.70 pb and backgrounds. Right: distribution of events with 0, 1 and > 2 6-tagged jets for the combined ^-|-jets channels, 
showing contributions from ti signal using a cross section of 7.93 pb and backgrounds. 

distribution for signal and background to the data. Systematic uncertainties are accounted for in the maximum 
likelihood fit by assigning a nuisance parameter to each independent systematic variation. Each nuisance parameter 
is modeled as a Gaussian probability density function with mean at zero and width corresponding to one standard 
deviation (SD) of the considered systematic uncertainty For a top quark mass of 172.5 GeV [i^, one obtains 
CTjt = 7.70j:o:7o (stat + syst + lumi) pb. 

The SM predicts that the top quark decays almost exclusively into a W boson and a 6-quark {t — )■ Wb). Hence 
besides using just kinematic information, the fraction of tt events in the selected sample can be enhanced using 
&-jet identification. One example to measure the tt cross section in a "counting" method is performed by the DO 
Collaboration analyzing 4.3 fb^^ of data. Final states with 3 and > 4 jets are explored. Each channel is further 
separated into events with 0, 1 and > 2 6-tagged jets, thereby obtaining 24 independent sets of data. Figure [1] (right) 
shows a comparison of the distributions in data with 0, 1 or > 2 6-tagged jets with the SM tt cross section, and the 
contributions from the different backgrounds. 

An iterative procedure is used to extract the tt cross section after 6-tagging. The fit of the tt cross section (utt) 
to data is performed using a maximum likelihood fit for the predicted number of events, which depends on att- To 
take into account each channel, the likelihood maximization procedure multiplies the Poisson probabilities for all 
channels j. The systematic uncertainties are incorporated in the fit using nuisance parameters [lO|, each represented 
by a Gaussian term. The result is a^j = 7.93t (stat + syst -I- lumi) pb. 

The CDF Collaboration, too, performs a "topological" measurement of the top pair production cross section in 
the £-|-jets channel with an integrated luminosity of 4.6 fb~^ using an artificial neural network (ANN) technique to 



discriminate between top pair production and background processes 12| . The ANN discriminant output is shown in 
Fig. m (left). 

The CDF Collaboration also performs a "counting" method analyzing the number of events with at least one 
6-tagged jet. The predicted number of events for each background process, along with the number of expected tt 
events at the measured cross section, is shown compared to data in Fig. [5] (right). To measure the tt cross section, 
a likelihood is formed from the data, the tt cross section, and the predicted background for that cross section. The 
result for the "topological" method is CTjj = 7.71±0.37(stat)±0.36(syst)±0.45(lumi)pb, the result for the "counting" 
method using 6-tagging is ctjj = 7.22±G.35(stat)±0.56(syst)±0.44(lumi)pb. The largest systematic uncertainties for 
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FIG. 2: Left: output of the ANN discriminant for £+> 3jets, showing contributions from tt signal using a cross section of 
7.71 pb and backgrounds. Right: distribution of events with > 1 6-tagged jets in the ^+jets channel, showing contributions 
from ti signal using a cross section of 7.22 pb and backgrounds. 



both experiments come from the measured luminosity and the 6-tag modeling in the simulation. Because 6-tagging 
is not used in the "topological" measurement, it is insensitive to the systematic uncertainty of the 5-jet tagging 
measurement. 

The dependence on the luminosity measurement and its associated large systematic uncertainty can be significantly 
reduced by computing the ratio of the ti to Z/7* boson production cross section, measured using the same triggers 
and dataset, and then multiply this ratio by the theoretical Z/7* production cross section [13|. In essence this 
replaces the luminosity uncertainty with the uncertainty on the theoretical cross section. The extracted 

tt cross section is a^j — 7.82 ± 0.38 (stat) ± 0.37 (syst) ± 0.15 (theory) pb for the "topological" method, a^j — 
7.32 ± 0.36 (stat) ± 0.59 (syst) ± 0.14 (theory) pb for the "counting" method with 6-tagging and (Tu = 7.70 ± 0.52 pb 
for the combination of both. The uncertainty of ±7% for the "topological" measurement is the most precise value 
determined by a single measurement to date. 

Figure [3] summarizes the measurements of the tt cross section in different decay channels performed by the CDF 
and DO Collaborations. All measurements agree with each other and agree with the SM predictions based on the full 
NLO matrix element 0| including soft-gluon resummation at next-to-leading logarithm (NLL) accuracy and at 
next-to- next-to-leading logarithm (NNLL) accuracy 0^ ■ The DO Collaboration has now analyzed all possible decay 
channels except final states with two hadronically decaying r leptons. In particular there is a new measurement in 
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FIG. 3: Summary of tt cross section measurements in different decay channels by the CDF (left) and the DO (right) Collabo- 
rations. The combination for each experiment is displayed, too. 
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FIG. 4: Predicted number of events as a function of the jet multiplicity (left) and Likelihood curve (right) for the measured 
tt+0 jet and tt+> 1 jet cross sections. 



the decay channel including a hadronically decaying r lepton and jets [ij] . The CDF Collaboration does not perform 
an explicit identification of t leptons but includes t final states in selections of missing transverse energy and jets [15] . 
This selection also recovers electrons and muons that were not selected in the i+jets and dilepton channels. It is 
also worth mentioning that both collaborations have achieved an uncertainty of 13% in measurements of the tt cross 



section in dilepton final states using a "topological" method [16| and a "counting" method using 6-tagging [17]. The 
CDF Collaboration combines the "topological" and "counting" measurements of the tt cross section in the ^+jets 
channel to the measurements in the dilepton fl^ and all-hadronic [3] channels and achieves an accuracy of 6%. The 
result is 

a^j = 7.50 ± 0.31 (stat) ± 0.34 (syst) ± 0.15 (theory + lumi residual) pb 
for a top quark mass of 172.5 GeV. 



B. Top Quark Pair + Jet Production 



An important test of perturbative QCD, as NLO effects play an important role in the calculation of the theoretical 
cross section, is the production of a pair associated with an additional hard jet (tt+jet). This process is interesting 
in its own right, because large fractions of the samples show additional jet activity and deviations from the SM could 
signal new physics such as top-quark compositeness. The first measurement of the cross section of tt+jet production 
has been performed with 4.1 fb^^ of collected data by the CDF Collaboration. The measurement is performed using 
^-tagged events in the ^+jets channel. A data-driven approach is used to predict the background content. The 
predicted number of events including the background, the tt+0 jets and the tt+> 1 jet contributions is compared to 
the data as a function of the jet multiplicity in Fig. |4] (left). A 2-dimensional likelihood is constructed from the data 
and prediction for events with three, four, or five jets to simultaneously measure the ti+jet and the tt without jet 
cross sections. This is shown in Fig. |4] (right). The measured result is dtj^jot — 1-6 ± 0.2 (stat) ± 0.5 (syst) pb which 
is in agreement with the SM prediction of (Tjj^jj,^ = 1.79j^Q'32 pb [l9| . 



C. Differential Top Quark Pair Production Cross Section 



Measurements of differential cross sections in the system test perturbative QCD for heavy-quark production, and 
can constrain potential physics beyond the SM. The transverse momentum (pt) of top quarks in tt events provides 
a unique window on heavy-quark production at large momentum scales. The DO Collaboration has performed a 
measurement of the ti cross section as a function of pr of top quarks for the first time [20]. Using a data set of Ifb^^ 
the £+iets final states with at least one 6-tagged jet are explored. 
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FIG. 5: Left: inclusive da/dpr for tt production (two entries per event for top and anti-top quark) in data (points) compared 
with expectations from NLO QCD (solid lines), from an approximate NNLO QCD calculation, and for several event generators 
(dashed and dot-dashed lines). The gray band represents the uncertainty on the QCD scale. Right: ratio of {\ / a)da / dpr for 
top quarks in tt production (two entries per event) to the expectation from NLO QCD. The uncertainty in the scale of QCD 
is displayed as the gray band. Also shown are ratios relative to NLO QCD for an approximate NNLO QCD calculation and 
of predictions for several event generators. In both plots inner and outer error bars represent statistical and total (statistical 
and systematic added in quadrature) uncertainties, respectively. 



Objects in the event are associated through a constrained kinematic fit to the — >■ WbWb — >■ ivbqq'b process. To 
obtain a background-subtracted data spectrum, the signal purity is fitted using signal and background contributions 
as a function of pr, and applied as a smooth multiplicative factor to the data. The reconstructed px spectrum is 
subsequently corrected for effects of finite experimental resolution, based on a regularized unfolding method [21]] 
using a migration matrix between the reconstructed and parton px derived from simulation. 

The measured differential cross section as a function of the pr of the top quark (using for each event the two 
measurements obtained from the leptonic and hadronic top quark decays), da/dpr , is shown in Fig. [5] (left) with the 
NLO QCD prediction 1221 . Also shown are results from an approximate next-to-NLO (NNLO) QCD calculation 23 [ 
and from the pythia [24| , the ALPGEN 



25| (matched with PYTHIA for parton showering |26l|) and the MC@nlo 27 



event generators. A shape comparison of the ratio of {1 / a)da / dpx relative to NLO QCD is shown in Fig. [S] (right). 
This shows that results from NLO and approximate NNLO QCD calculations and from the mc@nlo event generator 
agree with the normalization and shape of the measured cross section. Results from ALPGEN+pythia and pythia 
describe the shape of the data distribution, but not its normalization. 

The CDF Collaboration has measured the ti differential cross section in the £-t-jets channel with respect to the 
tt invariant mass, da/dMff, using an integrated luminosity of 2.7 fb~^ [3]. The tt invariant mass spectrum is a 
powerful test of higher order QCD calculations and also sensitive to a variety of exotic particles decaying into tt 
pairs. The tt invariant mass is reconstructed, using the four-vectors of the ^-tagged jet and the three remaining 
leading jets in the event, the lepton and the transverse components of the neutrino momentum, given by missing 
transverse energy. The expected contribution from the backgrounds is subtracted from the original M^j distribution. 
The resulting Mtt signal distribution suffers from resolution smearing and is corrected using a regularized unfolding 
technique, which also accounts for the longitudinal component of the neutrino momentum. The measured da/dM^f is 
shown in Fig. [5] (left). The result is consistent with the SM expectation, as modeled by PYTHIA with CTEQ5L parton 
distribution functions. However, the description by PYTHIA is not perfect. A significantly better description of the 
data is achieved in calculations where threshold logarithms up to NNLL order are resummed 0]. Fig. [S] (middle) 
shows a good agreement in the da/dMtt distribution with the same CDF data. 
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FIG. 6: Unfolded da/dMu distribution compared to pythia (left) and a NLO+NNLL QCD calculation (middle). Right: Mu 
distribution for SM ti production and backgrounds compared to data. 



III. SEARCHES FOR NEW PHYSICS IN TOP QUARK PAIR PRODUCTION 

The fact that the top quark has by far the largest mass of all known elementary particles suggests that it may 
play a special role in the dynamics of electroweak symmetry breaking. This makes searches for new physics in the 
top quark sector very attractive. In the following one recent example for a search for new physics in the production 
of top quarks is presented: a search for ti resonances which explores the invariant ti mass distribution, too. 

One of the various models incorporating the possibility of a special role of the top quark in the dynamics of 
electroweak symmetry breaking is topcolor [29| , where the large top quark mass can be generated through a dynamical 
ti condensate, X, which is formed by a new strong gauge force preferentially coupled to the third generation of 
fermions. In one particular model, topcolor-assisted technicolor 30], X couples weakly and symmetrically to the first 
and second generations and strongly to the third generation of quarks, and has no couplings to leptons, resulting in 
a predicted cross section for ti production larger than SM prediction. 

The CDF and DO collaborations presented updated model-independent searches for a narrow- width heavy resonance 
X decaying into ti using 2.8 fb~^ of data in the all-hadronic final state (CDF) 31] and 3.6 fb~^ of data analyzing 
£-|-jets final states (D0) [3^. A new search for resonant production of ti is performed by the CDF collaboration 
in the .^-l-jets channel using a data set of 4.8 fb^^. The search for resonant production is performed by examining 
the reconstructed ti invariant mass distribution which is shown in Fig. [5] (right). A matrix element reconstruction 
technique is used and for each event a probability density function (pdf) of the ti invariant mass is sampled. These pdfs 
are used to construct a likelihood function, whereby the cross section for resonant ti production, given a hypothetical 
resonance mass and width, is estimated. One compares the 95% C.L. points for these posterior probabilities for the 
data with those expected from SM sources only and determines how likely it is for the SM to fluctuate to the data for 
each mass point. A very unlikely fluctuation could then indicate the presence of new physics in the sample. There is 
no evidence of resonant production of ti candidate events. Therefore within a top-color-assisted technicolor model, 
the existence of a leptophobic Z' boson with Mz' < 900 GeV and width around Tz' ~ Q.012Mz' is excluded at 95% 
CL. 



IV. SINGLE TOP QUARK PRODUCTION 



A. Single Top Quark Production Cross Section and Vtt 



Single top quark production serves as a probe of the Wtb interaction [33[ , and its production cross section provides 



a direct measurement of the magnitude of the quark mixing matrix element Vtb without assuming three quark 



generations 3J]. However, measuring the yield of single top quarks is difficult because of the small production rate 



and large backgrounds. At the Tevatron single top quarks are produced by either a i-channel exchange of a virtual 
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FIG. 7: Output discriminants for single top quark observation. Displayed are the CDF super-discriminant output (upper left), 
the CDF discriminant output for the analysis using a missing transverse momentum and jets selection (upper middle) and the 
DO BNN discriminant (lower left). A summary of all existing measurements and the Tevatron combination is shown (lower 
right). 



W boson which combines with a highly energetic h quark to produce a top quark, or by an s-channel exchange of a 
far off-shell W boson which decays to produce a top quark and a b antiquark. 

In March 2009, the CDF and DO Collaborations reported observation of the electroweak production of single top 



quarks in pp colhsions at ^ = 1.96 TeV based on 3.2 fb"^ (CDF) [35|, Isg and 2.3 fb^^ (D0) of data. Both 
Collaborations used events containing an isolated electron or muon and missing transverse energy, together with jets 
where one or two of the jets were required to originate from the fragmentation of b quarks. The data are devided 
into independent sets (subchannels) depending on lepton flavor, jet multiplicity, number of 6-tags and in the case of 
DO the run period. 

CDF and DO each combine many variables using different multivariate analysis methods such as Boosted Decision 
Trees, Neural Networks, Bayesian Neural Networks, Matrix Elements and Likelihood functions to increase the sepa- 
ration power between signal and background. The discriminant outputs of each multivariate analysis are combined to 
one discriminant taking the correlations into account. CDF combines 8 subchannels into a super-discriminant using 
a neural network trained with neuroevolution 13811 which is shown in Fig. [7] (upper left). The CDF Collaboration has 
an additional independent search channel 35|, |39| that is designed to select events with missing transverse energy 
and jets and is orthogonal to the other channels. It accepts events in which the W boson decays into hadronically 
decaying t leptons and recovers events lost due to electron or muon identification inefficiencies. The discriminant 
output of this analysis - shown in Fig. [7] (upper middle) - is combined to the super-discriminant to obtain the final 
result. DO combines 12 subchannels into one BNN discriminant displayed in Fig. [7] (lower left). 

CDF measures a cross section of cr(pp tb + tqb + X) = 2.31q'5 assuming a top quark mass of 175 GeV 
and DO measures a cross section of a{pp — > tb + X, tqb + X) = 3.94 ± 0.88 pb at a top quark mass of 170 GeV. Both 
measurements correspond to a 5.0 standard deviation (SD) significance for the observation. They are in agreement 
with the SM predictions. Both results are translated into a direct measurement of the amplitude of the CKM matrix 
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element Vtb without making assumptions on the number of quark generations and the matrix unitarity. CDF obtains 
\Vtb\ = 0.91 ±0.13, DO derives \Vtb\ = 1.07 ±0.12. 
AU measurements mentioned above 



35l 1371 l39j and a measurement by the DO CoUaboration analyzing the r+jets 



41 



final state by explicitly identifying hadronically decaying r leptons [40] , agree with the SM predictions 
can be seen in Fig. [7] (lower right). The Tevatron combination (still excluding the r+jets final state) |4 
single top cross section of 



4^. This 



gives a 



a{pp-^ tb + X, tqb + X) 
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pb 



and a CKM matrix element of 



\Vtb\ = 0.88 ±0.07. 



B. s-Channel and t-Channel Single Top Quark Production 

In the observation analyses the combined s + t channel single top quark cross section was measured, assuming the 
SM ratio of the two production modes. This ratio is modified in several new physics scenarios, for example in models 
with additional quark generations, new heavy bosons, fiavor-changing neutral currents, or anomalous top quark 
couplings. Therefore it is interesting to remove this constraint and to use the t-channel characteristics to measure 
the t-channel and s-channel cross sections simultaneously which provides a t-channel measurement independent of 
the s-channel cross section model. 

The DO Collaboration reported direct evidence for the electroweak production of single top quarks through the 



t-channel exchange of a virtual W boson alone 44|. The measured cross section is a{t — channel) = 3.14to'gi pb, has 
a significance of 4.8 SD and is consistent with the SM prediction. This is the first analysis to isolate an individual 
single top quark production channel. Both the CDF and the DO Collaborations extracted the s-channel (cts) and 
t-channel {at) cross sections simultaneously relaxing the SM ratio assumption for Os/dt 3^ 44 1. The results agree 
with the SM predictions. 



C. Searches for New Physics in Single Top Quark Production 

The single top quark production channel offers a large variety for sensitive searches for new physics beyond the 



SM. One example are searches for fiavor changing neutral current (FCNC) couplings [45|, |46|. Recently, the DO 
Collaboration searched for single top quark production via FCNC top^luon-quark couplings Ktgu and Ktgc in a 
sample corresponding to 2.3 fb~^ of integrated luminosity collected [46|. Events containing a single top quark 
candidate with an additional jet are selected. Separation between signal and background is obtained using Bayesian 
neural networks (BNN) which are trained for each choice of lepton flavor (electron or muon), jet multiplicity (2, 3, 
or 4), and data-taking period separately, twelve in total. Fig. [8] (left) shows the comparison between background and 
data for all twelve BNN discriminants combined. Since the data are consistent with the background expectation, 
upper limits on the FCNC cross sections and couplings are set using a Bayesian approach [47]. Following the analysis 
strategy of previous work [4^ , a two-dimensional Bayesian posterior density for the cross sections and for the square 
of the FCNC couplings is formed, using the BNN distributions for data, background, and signals. Fig. |S] (right) 
shows the posterior density as a function of the FCNC cross sections Utgu and utgc- Since there is agreement with 
the SM, limits are set on the couplings of Ktgu/ ^ < 0.013 TeV~^ and Ktgd !^ < 0.057 TeV^^, A being the scale of the 
new interactions which generate these couplings (of order 1 TeV), without making assumptions about the tgc and 
tgu couplings, respectively. 
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FIG. 8: Left: comparison of the background model to data for the FCNC discriminant summed over all analysis channels. 
The bins have been ordered by their signal to background ratio and the FCNC signals are each normalized to a cross section 
of 5 pb. Right: Bayesian posterior probability as a function of the atgu and atgc cross sections. 



V. CONCLUSIONS 

Recent highlights in top quark production analyses from the Tevatron coUider are reviewed. Among many impres- 
sive results, the observation of single top quark production via the electroweak interaction and the direct measurement 
of the CKM matrix element Vtb by the CDF and DO Collaborations is outstanding. 

Analyzing datasets corresponding to an integrated luminosity of up to 5.7 fb^^ the Tevatron experiments can 
perform high precision measurements. As an example the top quark pair production cross section is extracted with 
an accuracy of 6%. The high statistics in ti samples allows to analyze differential cross sections for the first time 
with a high accuracy resulting in powerful tests of higher order QCD calculations. In general the inclusion of soft 
gluon resummations at the NNLL level into NLO QCD calculations improves the description of the data. Both top 
quark pair production and single top quark production channels allow sensitive searches for new phenomena beyond 
the SM. So far, there is no hint for new physics beyond the SM in the top sector. Excellent prospects for top quark 
physics investigations can be expected with more integrated luminosity at the Tevatron. 
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